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Abstract 

We consider shock measures in a class of conserving stochastic parti- 
cle systems on Z. These shock measures have a product structure with a 
step-like density profile and include a second class particle at the shock 
position. We show for the asymmetric simple exclusion process, for the 
exponential bricklayers' process, and for a generalized zero range process, 
that under certain conditions these shocks, and therefore the second class 
particles, perform a simple random walk. Some previous results, including 
random walks of product shock measures and stationary shock measures 
seen from a second class particle, are direct consequences of our more gen- 
eral theorem. Multiple shocks can also be handled easily in this frame- 
work. Similar shock structure is also found in a nonconserving model, 
the branching coalescing random walk, where the role of the second class 
particle is played by the rightmost (or leftmost) particle. 

Keywords: Interacting particle systems, second class particle, shock measure, 
exact solution, asymmetric simple exclusion, zero range process, bricklayers pro- 
cess, branching coalescing random walks 
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1 Introduction 

On a macroscopic level driven diffusive systems are often described by a set 
of conservation laws for the densities. These hydrodynamical equations are in 
general nonlinear PDEs, which can develop singularities in the solution. Shocks 
are discontinuities in these weak solutions, which travel with a speed also known 
as the Rankine-Hugoniot velocity. Whereas the large scale continuous descrip- 
tion of shocks is well established [T] , much less is known about the microscopic 
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structure and dynamics, which has become a subject of intense investigation in 
recent years [3 11 [3 11 [3 H 1 [13 [III HI III III] • 

In order to make the first steps in this direction one has to define the position 
of shocks on the lattice scale, which is already a nontrivial task in general. It is 
well known that second class particles, which move stochastically and follow the 
trajectories of density fluctuations, are attracted by shocks and therefore serve 
as good markers for the shock position. Derrida et al. in [5] derive the time 
invariant shock measure in the asymmetric simple exclusion process (ASEP) as 
seen from such a second class particle. They observe that when a certain condi- 
tion holds for the asymmetry and the limiting densities, the invariant measure 
becomes a Bernoulli product measure with a simple step-like density profile. 
Similar product shock structure, as seen from a second class particle, was found 
later in another stochastic lattice model, the exponential bricklayers' process 
(BLP) iJ. 

A different approach, initiated by Belitsky and Schiitz, attempted to capture 
not only the structure but also the microscopic dynamics of shocks. In [6] 
they show that under the same condition as in [5] there is a family of product 
measures /x^, fc € Z with a step-like density which evolve into linear combinations 
of similar measures, and the interpretation is that the shock position k performs 
a simple random walk. The random walking shocks were shown to exist later in 
the exponential BLP [7] too. The advantage of this description is that it doesn't 
use second class particles therefore it can be applied also in cases where second 
class particles cannot be defined (or the number of them is not conserved). An 
example of such case is the branching coalescing random walk. Although this is 
a non- conserving system, there are shocks with similar structure and evolution 
here as well [8]. It is interesting that random walking shocks have also been 
found in systems with more than one conserved quantities 0. 

These two types of results naturally raise the question whether the second 
class particle itself, attracted by the shock, performs a simple random walk. In 
this paper we give an answer to this question by considering shock measures with 
second class particles at the shock position. The idea of considering such shock 
measures appeared in the context of the ASEP with open boundaries in [8], 
where a conjecture is formulated saying that the above random walk property 
should hold for the shock measures with second class particles too. In a fairly 
general framework we show the random walk property for these measures in the 
ASEP, in the exponential BLP and in a generalized zero range process (GZRP) 
(where negative particle numbers can also occur). While our result clearly 
shows the simple random walk of second class particles, the results of [2l|5l|6j[7] 
for product shock measures also follow. Notice that the existing stationary 
product distribution results did not include any random walk dynamics, and 
the random walk results did not include the second class particle. Hence our 
result is genuinely new, and also connects the two types of arguments. 

The diffusion coefficient of a general shock in ASEP was computed and the 
diffusive behavior was investigated by Ferrari and Pontes [3]. They used the 
second class particle as a characterizing object for the shock location, and their 
result on how the shock location depends on the initial configuration of ASEP 
made it possible to generalize the diffusive scale-results to the case of multiple 
shocks with second class particles in Ferrari, Pontes and Vares [4]. 

Multiple shocks, i.e., several steps in the density profile have also been stud- 
ied in [5] and [J. In this case the exact microscopic description is given, and 



2 



involves several shock positions, which perform an interacting simple random 
walk. Due to the attraction of these "micro-shocks" they form a bound state 
with a finite width and can be considered as a single shock with a more com- 
plex structure. It is interesting that while the result for multiple shocks is a 
direct generalization of that for a single shock in the BLP [7] , in the ASEP such 
a naive generalization does not hold. In order to be able to handle multiple 
shocks in the ASEP, extra particles or vacancies had to be introduced at the 
shock position fBi. Our description with the second class particles at the shock 
positions explains the interaction between shocks in a very natural way, without 
the need of artificial particles. In fact we show that previously known forms of 
random walking shocks can be obtained as an appropriate mixture of the two 
marginals of our coupling shock measure. 

It is important to note that the existence of random walking shocks is closely 
related to the exact solvability of these particle systems in an open geometry. 
The matrix product ansatz is a method widely used for finding stationary states 
of stochastic lattice models [15] . Recently it has been shown jlOj [E] that the 
occurrence of a single random walking shock implies the existence of a two- 
dimensional representation of the quadratic algebra appearing in the matrix 
product ansatz. Similarly, multiple shocks correspond to other finite dimen- 
sional representations [TnillUlIS]- Moreover, as a very recent progress, utilizing 
the single-particle properties of random walking shocks, the Bethe ansatz has 
been successfully applied [13] in the open-boundary ASEP for the evaluation of 
the spectrum and also for the current large deviation function. 

The paper is organized as follows. In section[5]we introduce a wide family of 
stochastic particle systems on Z with nearest neighbour jumps, and summarize 
their basic properties, such as stationary product measures and hydrodynamic 
limit. A few specific examples are considered in more detail. A surface growth 
interpretation can also be given to these models which is sometimes more nat- 
ural. In section [3] we formulate our main results for random walking single and 
multiple shocks, proofs are given in section UJ A different model, the branching 
coalescing random walk is investigated in section [S] This model is not in the 
family considered before. We prove a statement here which is similar in spirit 
to that of section [3J The second class particle is replaced here by the rightmost 
(or leftmost) particle. 

2 A family of models 

The class of stochastic interacting systems we consider here appeared several 
times in the literature, we repeat a description recently formulated in |16| . The 
class is a generalization of the so-called misanthrope process. We use a surface 
growth interpretation, but many members of this class can be understood in 
terms of particles jumping on the one dimensional lattice. For — oo < w™™ < 
and 1 < ijj™-'^'^ < oo (possibly infinite valued) integers, we define the single-site 
state space 

/ : = {z e Z : w""" -l<z< uj"'"'' + 1} 
and the configuration space 
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Figure 1: Shown is the relation between the surface (above) and the particle 
(below) interpretation. The empty circle denotes an antiparticle. 



For each pair of neighboring sites i and i + 1 of Z, we consider a column 
built of bricks above the edge (i, i + 1). The height of this column is denoted 
hy hi. A state configuration co G fl has components uji = /i^-i — hi S /, being 
the negative discrete gradients of the height of the "wall" . For co^il and i ^ j 
let w''-' be the configuration with components 



(2.1) 



Also, define, for a vector h of heights, /i*^ and h^^ by 

1 ../o, for k i, , I f hk, for k ^ i, 

/ii, = < , ^ , h,* 





1, for k = i, ^ /ifc — 1, for k = i. 

The continuous time evolution is described by jump processes. A brick can 
be added: 

> with rate p{uji, w^+i), 



LO > LJ 

h h'^ 



or removed: 



h — > h'^ 



with rate q{uJi, i-^i+i)- 



Conditionally on io(t), these moves are independent. See Figure [T] for visualiza- 
tion. We impose the following assumptions on the rates: 

• The rates must satisfy 

^(0;""", cj"^^'^) = qiuj"'^'', ■ ) = g( • , ^""") = 
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p and q are non-zero in all other cases, or one of them is the identically 
zero function (totally asymmetric case). 

• In order to provide a smoothening effect in the dynamics we assume mono- 
tonicity in the following way: 

^2 2) P{z + l,y)>p{z,y), p{y,z + l)<p{y,z) 

q{z + 1, y) < q{z, y), q{y, z + l)>q{y, z) 

for y,z,z + lEl. This property has the natural interpretation that the 
higher neighbors a column has, the faster it grows and the slower it gets 
a brick removed. Our model is hence attractive. 

• We are going to use the product property of the model's translation- 
invariant stationary measure. For this reason, similarly to Cocozza-Thi- 
vent [T71 . we need two assumptions: 

— For any x, y, z £ I 

p{x, y)+p{y, z)+p{z, x) 
(2.3) +q(x, y) + q{y, z) + q{z, x) ^p{x, z)+p{z, y)+p{y, x) 

+ q{x, z) + q{z, y) + q{y, x). 

— There are symmetric functions Sp and Sq, and a common function /, 
such that /(w™'") = whenever w™'" is finite, and for any y, z G / 
(2.4) 

p(y, z) = Sp(y, z + 1) • /(y) and q{y,z)^Sg{y+l,z)-f{z). 
Condition (12. 2p implies that / is non-decreasing on /. 

• In order to properly construct the dynamics, restrictive growth conditions 
might be necessary on the rates p and q in case of an unbounded single-site 
state space /. We comment on this below. 

At time the interface mentioned above is described by w(t) . Let (ys : — ;> R 
be a finite cylinder function, i.e., depends on a finite number of Ui values. 
The growth of this interface is a Markov process, with the formal infinitesimal 
generator L: 

(L(^)(w) ^p(w,, w,+i) • [(^(0;*''+^) - ^{ui)] 

(2.5) 

The construction of dynamics is available in the following situations. Several 
models with bounded rates are well understood and can be handled via the Hille- 
Yosida Theorem, see Liggett p[l] . When the rates p and q grow at most linearly 
fast as functions of the local w values, then methods initiated by Liggett and 
Andjel lead to the construction of some zero range type systems (Andjel [19], 
Liggett [5^, Booth and Quant [2I1I22])- The totally asymmetric zero range and 
bricklayers' processes with at most exponentially growing rates are constructed 
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in Balazs, Rassoul-Agha, Seppalainen and Sethuraman 123| . See the definition 
of zero range and bricklayers' processes below. 

We assume that the existence of dynamics can be established on a set of 
tempered configurations 17 (i.e. configurations obeying some restrictive growth 
conditions), and we have the usual properties of the semigroup and the generator 
acting on nice functions on this set. We also assume that fl is of full measure 
w.r.t. the stationary measures defined in Section [521 Questions of existence of 
dynamics are not considered in the present paper. 

2.1 Examples 

Many well-known nearest neighbor processes belong to this class, see e.g. [TB] 
for a more complete treatment. Here we only list those which we consider later 
in more detail. 

• The asymmetric simple exclusion process (ASEP) introduced by 
F. Spitzer ^ is characterised by w""" = 0, cj'"'''' = 1, /(z) = l{z = 1}, 

Sp{y, z) =p-l{y = z = 1} and Sg{y, z) ^ q ■ l{y = z ^ I}, 

where p > q are non- negative reals adding up to 1 (see (12. 4p ). In this case 

P{y, z)=p-l{y=l, z = 0} and g(y, z) = g • l{y = 0, z = 1}. 

Here uji G {0, 1} is the occupation number for site i, p{LOi, i^i+i) is the 
rate for a particle to jump from site i to i + 1, and q{uJi^ ^i+i) is the rate 
for a particle to jump from site i + 1 to i. These rates have values p and g, 
respectively, whenever there is a particle to perform the above jumps, and 
there is no particle on the terminal site of the jumps. Conditions (12. 2p 
and (j2.3|) are also satisfied by these rates. 

• Totally asymmetric zero range processes are included by an arbi- 
trary nondecreasing function / : Z — >■ R+ , 



In its original form, the totally asymmetric zero range process is a particle 
system with oji particles at site i, a particle jumps from z to i -j- 1 with 
rate /(w,). The setting w""" = 0, w™'*'' = oo, and so /(O) = would 
correspond to this situation. In the sequel it will be important for us 
to allow for negative values of u as well, which comes naturally in the 
surface representation. We shall refer to this class of models with a;™'" = 
— oo, ijj™-^^ — oo, and f{y) > for all y S Z, as generalized zero range 
processes. Conditions (12.21) and (j2.3p trivially hold for the rates. 

— As a special case, the generalized totally asymmetric exponen- 
tial zero range process (we will simply refer to it as GZRP) is 

obtained by p{y, z) — /(y) — qP'^v^'^/'^) with a /? > parameter. 
Omitting the constant — /3/2 would simply correspond to a change 
of timescale but would bring in some unwanted factors in our final 
result. 



Sp(y, z) = 1 
P{V, z) = f{y) 



and 



and 



Sq{y, z) = 0, 
q{y, z) = 0. 
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• Totally asymmetric bricklayers models. Let w""" — —oo, cj™^^ — cxd, 
/ : Z — > M+ non-decreasing, also having the property 

/(z) • /(I - z) = 1 for ah zeZ. 

The values of / for positive z's thus determine the values for non-positive 
z's. Set ^ 

Sp{y, z) = 1 + jj-jj^, Sq{y,z) = 0, 

which results in 

Piy, z) = f{y) + f{-z), q{y,z) = 0. 

This process was first represented by bricklayers standing at each site i, 
laying a brick on the column on their left with rate /(— Wi) and laying a 
brick to their right with rate f{uJi), hence the name. Conditions (j2.2p and 
(j2.3p hold for the rates. 

— As a special case, the totally asymmetric exponential brick- 
layers process (we will just abbreviate as BLP) is obtained by 
/(z) = Qpi^-^/^) with a /3 > parameter. 

Note that the BLP is a symmetrized version of the GZRP in the following 
sense. For a GZRP with a given l3 one can consider its counterpart by a 
space reflection in the surface representation. The symmetric combination 
of these two processes (which is obtained by taking the sum of the two 
generators) gives the corresponding BLP (in the particle representation a 
particle-antiparticle transformation should follow the space reflection). 



2.2 Translation invariant stationary product distributions 

We now present some translation invariant stationary distributions for these 
processes. For many cases it has been proved that these are the only extremal 
translation-invariant stationary distributions. Following some ideas in Cocozza- 
Thivent [ITj, we first consider the non-decreasing function / of (|2.4p . For / 3 
z > we define 

z 

1 

II f{y) 

y=z+l 

finally /(O)! : = 1. Then we have 

/(z)!./(z + l) = /(z + l)! 

for all z e /. Let 

f lim log(/(z)) , if w'"'''' = oo 
I oo , else 



while for / 9 z < let 
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and 

f lim log(/(-z)) , if w""'" = -oo 
~ 1 —03 , else. 

By monotonicity of /, we have 9 > 9. We assume 9 > 9. With a generic real 
parameter 9 £ {9., 9), which is often referred to as the chemical potential, we 
define the partition sum as 

Let the product-distribution /i^ have marginals 

1 e« 



if z e /, 



(2.6) ^Ji'{z)=^i^'{^■.Lo, = z}■.= {Z{9) f{z)l 

if z ^ /. 

Then the product distribution fi^ is stationary for the process generated by 
We define the function 

(2.7) p(0):=^log(Z(0))-^z./(z), 9<9<9. 

zei 

This is the density of particles, and it is quite easy to see that this is a strictly 
increasing function. Its inverse will be denoted by 9(g). Due to the bijection it 
is always possible to choose a suitable characterization of the stationary state 
either with the chemical potential 9 or with the corresponding density g. 

As for our examples, the ASEP has a Bernoulli product stationary distribu- 
tion of which the density parameter is 

(2-8) PiO) = 

In the case of both our exponential GZRP and BLP examples, computing 
the factorials in (12.61) results in the discrete Gaussian 



from which 



00 / \ 2 



Z——00 



While no explicit form is available for the partition sum Z{9), the identity 

(2.9) Zi9-I3)^e^/^-^ ■Z{9) 
can easily be shown. Via ^2.7\i this implies 

(2.10) pi9 - (3) ^ pi9) ~ 1 

for the exponential GZRP and BLP processes. Notice also that 

(2.11) A^(y) = /(2/-t-l). 
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2.3 Hydrodynamics, very briefly 

In preparation for stating the result, we first briefly mention that it is beheved, 
and in many cases proved, that models of our family satisfy a conservation law 
of the form 

(2.12) dTp{T, X) + dxn{p{T, X)) = 

in the Eulerian scaling^ with the density p being a function of the rescaled time 
and space variables T and X, and Tiig), the flux function, being the expected 
net current in the stationary measure with density g (recall (12.71) '): 

nig) - WJ^' ^) - 9(2/' ^)] • m'^^^(2/)/^^'(^). 

See e.g. Rezakhanlou [33] or Bahadoran, Guiol, Ravishankar and Saada [33] for 
details. 

At time zero let the density proflle be a step function: p{0, X) = g for 
X < and p{0, X) = X for X > 0. This initial condition corresponds to a 
shock if g < X and Ti is concave, or if > A and f-L is convex. The entropy 
solutions of (|2.12l) in these cases are rigid translations of the shock with the 
Rankine-Hugoniot velocity 

(2.13) ^^mX) njg)^ 

X - g 

From this formula one can check that multiple shocks eventually meet and merge 
into a single shock having the leftmost and the rightmost of the initial density 
values on its left and right sides, respectively. 
It is easy to see that for the ASEP we have 

'Hig) ^{p-qj- g{l - g) 

which is a concave function (p > q). It is also straightforward from the deflni- 
tions that 

n{g) - e^^e) and U{g) = e^^^^ + e^^^^' 

for the GZRP and BLP examples, respectively. Due to the lack of an explicit 
formula for d{g), it is a nontrivial fact that both these flux functions are convex 
when the rate functions / are convex [27] . 



2.4 The second class particle 

Attractivity (|2.2p makes it possible to define the second class particle. Let 5^ be 
the vector of components 1 for site j, and for all other sites in Z. Let a; S 51 
such that ujj < w™'*'^, and 

(2.14) C: = w+^^'. 

We say that there is a second class particle at site j. We let the pair (w, Q 
evolve in the basic coupling. That is, assuming and making use of (j2.14p . the 
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effect of tfie coupled generator (which, we also denote by L) for a finite cylinder 
function (/? of a pair (w, ^) is 

(2.15) 

i—a— 1 

+ E p(^„ . [^(^^•'+\ - '/'(^^, C)] 

+ p(c.,_i, c., + 1) . [ifiuj^-''^, C'"-^'-'") - ^i^, 0] 

(2.16) + [p(w,-i, Lo,) - p(c^j_i, c^j + 1)] ■ [ifiui^-''', C) - ^i^, 0] 
+ p(w„ c.,+1) . [^{ui=^^+\ e'^') - ifiu;, 0] 

(2.17) + [piuj, + 1, ^,+i) - p(a.j, ^,+i)] • [^{ui, - ^{ui, 0] 

i—a—1 
b 

+ q{u,^u uJj) ■ [^(^"-\ C""') - ^i^, 0] 

(2.18) + [q{uj,-u + 1) - q{uj,^i, u,)] ■ [^{^, C""') ~ fi^, 0] 
+ q{u;, + 1, uj,+,) ■ [^(^■'"+l^^ f +1'^) - ^{ui, 0] 

(2.19) + [q{u,, UJ,+i) - q{uj, + 1, Wj+i)] • [ip{u£+^-\ - 0] , 

where Lp depends on the configuration over the sites a . . .b. We assume that 
a<— 1<1<6. This generator gives the correct marginal evolution for each 
Lo and C- Only steps p.l6p - p.l9[) influence the second class particle. In fact 
these steps result in jumps of this particle, hence the single second class particle 
is conserved for all times. 

Constructing the dynamics of a coupled pair with any number of second class 
particles can be done along the same lines, and their number is again conserved, 
see the generator for that case in Section [3?2] 



3 Shock measures with second class particles 

We now define product shock measures with a single second class particle on a 
coupled pair of processes. Later on multiple shocks and second class particles 
will be considered, but the case of one second class particle is much simpler, so 
we first demonstrate it here. The marginals will be like the (|2.6p stationary ones, 
except for one site which has the second class particle. Define, for ^ < < ^, 
the one-site marginal on I x I hy 



0, if y 7^ z, 
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where is the stationary marginal (|2.6p . Let also /t be a measure on / such 
that /i(a;™^'^) = if w'"^'' is finite. Define 




0, otherwise. 



With these marginals and with 9 < 9, a < 9 we define the product measure 
(3.1) — (S)^'^'^^^" 

of marginals on the left of site j, v at j, and v'^ on the right of site j. 
Therefore, v_j is a measure on coupled pairs with exactly one second class particle 
at site j. 

In the case of the ASEP it is more convenient to use the densities as param- 
eters instead of fugacities. Here, with a slight abuse of notation, the results are 
expressed in terms of p = p{9) and A = p{a) being the left and right densities 
correspondingly, according to (j2.8p . 



3.1 Results for a single shock 

The main result is on the time-evolution of some particular distributions of a 
coupled pair C(0) of some particular models. We use the semigroup nota- 

tion for the evolution of such a distribution, so vS(t) denotes the distribution of 
the pair at time t if the initial distribution at time was i/. Recall the definition 

Theorem 3.1. The identity 

(3.2) -^^v,^S{t) _^^P- - + Q ■ - ^^.] 



t=o 



holds, in the sense of test functions, in the following special cases with the 
following parameters P and Q: 

• For the ASEP, if the relation 

(3.3) ^ = ^ 

e(l-A) q 

holds between the densities and the asymmetry, and 

(3.4) A(0) = 1. 



In this case 



(3.5) 



P = ■ p = — ■ q = {\ — X)p + Xq and 

1- g g 

Q = Y~rj ■ Q = J ■ p ^ - 9)1 + QP- 



11 



• For the exponential GZRP and BLP as defined in Section \2. 11 if the rela- 
tion 



(3.6) 



0-(T = l3 



holds between the parameters, and 



(3.7) 



A(2/)=m"(2/)- 



In this case we have 



P = - 



and 



Q = 
Q = e 



for the GZRP, 
for the BLP. 



(3.8) 



P = - e"" 



and 



— a 
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— e 



The only way to have a second class particle at j in ASEP is Uj = and 
Q = 1. This fact is of course reflected in the form of j>. This is not the case for 
GZRP and BLP, where the marginals for ujj and Q are not a priori restricted 
in any way. The result turns out to require these marginals to be /i*^ and /i^, 
respectively (this latter being a consequence of (|2.1ip ). 

Remark 3.2. As in IQj and '7j, (j3.2[) has a natural random walk interpretation: 
the shock measure v_j , with the second class particle in the middle, performs an 
asymmetric simple random walk. The drift P ^ Q of this walk agrees with the 
Rankine Hugoniot velocity (j2.13p . This can be seen for the ASEP from the 
definitions and using the relation (|3.3p . For the exponential GZRP and BLP 
models the relation (|3.6p together with (|2.10p can be used to check this. 

We remark here that a theorem, similar in spirit, will be shown in Section 
15.21 for the branching coalescing random walk model. 

In 6 a shock measure with density g for sites at or left of j and A on 
the right of site j was first considered. It was shown there that this structure 
performs a random walk with the above jump rates. In view of our result, 
mixing the two marginals of leads to [6j's shock measure: flip a biased coin 
initially, and with probability A follow the distribution of the upper process 
C while with probability 1 — A follow the lower process w. Later on, in their 
Theorem 2, Belitsky and Schiitz in [B] put particles at the positions of the shocks 
to describe interaction (see our Section [X^ for interaction of shocks), and they 
also remark that these are not the second class particles. However, their setting 
is also obtained from our results by considering the second marginal Q of v_y 
Therefore, their particles at the shock positions can also be viewed as second 
class particles. Later on, the idea of random walking shocks with second class 
particles at the shock positions appeared in the form of a conjecture in |8j in 
the context of the ASEP with open boundaries. 

We recover the results of 7' if we only consider the first marginal a; of v_j 
for the exponential BLP. 

3.2 Results for multiple shocks 

While we tried to keep notations as simple as possible for a single shock, multiple 
shocks seem to require a more complicated treatment. We describe below and 
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use a formalism capable of handling multiple shocks with second class particles. 
First we rewrite (|2.15|) in the general case, only assuming uji < Q for each i € Z: 

(3.9) {L^){^, C) = E C.+i) • C"+') - ^(^, 0] 

i—a — l 

(3.10) +[p{0J,, LJ,+,) - p{lO,, C^+l)] ■ C) - ¥'(^, 0] 

(3.11) +[p{a, Q+i)-p{io,, c+i)] • [^(^, C''+') - ^(^, 0] 

(3.12) +[g(c., c,+i) - 9(0, ^.+i)] • H^^, c^'n - ^(^, 0] } 

(3.13) +[<z(c^„ w,+i) - q{Q, w.+i)] • b(^'+''\ - ^{^, C)]- 

Again, the above differences of jump rates are nonnegative, and marginally both 
a; and C evolve according to the original dynamics. Moreover, the generator 
keeps LOi < Q for each site. 

Let m > be an integer and 6 < a < 9 a parameter value. Define the 
one-site marginal of a coupled pair w, C by 



A"^'"(y), ifz = y + m, 
0, ii z =/= y + m. 



In our examples this marginal will coincide with the stationary marginal (|2.6p 
if m = 0. When m > 0, it describes a distribution on a site with m second class 
particles; uj ~ ( ~ m then has distribution /i'^''". We require that /i"^'™ gives 
probability zero on values y > co'^'^^ — m when ui'^^^ is finite. Now with a vector 
m of nonnegative integer components rrii and a parameter vector ct, define the 
product measure 

(3.14) : = (g)J/'^•.'"• 

^ez 

on pairs (w, (). This measure describes m, second class particles at site i with 
marginal for Wi = ^ — and, in our examples, stationary marginals 

^(Ti.o _ ^<7i other sites with no second class particles. We assume a finite 

number of second class particles in the system that is, ^ < oo. 

iez 

In the ASEP case we still use (|2.8p and switch to the Qi parametrization. 
We also adopt the notation (|2.ip for m*'^, while ct*'-' will be some modified 
parameter vector which we specify below in the theorem: 

Theorem 3.3. The identity 

d 



(3.15) 



t=o 

-|-(5(TOi, TTT-i+l, ai, ffi+l) ■ [v- 



holds, in the sense of test functions, in the following special cases with the follow- 
ing parameters a'''^'^^ , q^^^-^ , P{mi, rui^i, ai, Ci+i) and Q{mi, rui+i, Oi, <Ji+i): 
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For the ASEP, if mi =0 or 1, and the relation 

n CI n\ f 1' j/"lj+l=0, 

3.16 — ^ I y 

holds between the densities and the asymmetry, and 



(3.17) /i^'™(0) = 1 - ^"^'"(1) = 

In this case 



1- g, ifm = 0, 
1 , if m — 1. 



(3.18) ,r^=r- 

[Qi-i, forj = i, 

(3.19) ,^'^ = 1^- f'^^^^' ^ 

( ft+i, /or J = I, and 

(3.20) P(mi, m^+i, gi, gi+i) = mj(l - m^+i) . [(1 ft+i)p + ft+i?], 

(3.21) Q(TOi, m^+i, ft, ft+i) = (1 - mi)mi+i ■ [(1 - ft)g + ftp]. 

• For the exponential GZRP and BLP as defined in Section \2.1[ if the rela- 
tion 

(3.22) cFi^i - ai = I3mi 
holds between the parameters, and 

the stationary marginal (|2.6p . regardless of m. In this case we have 



(3.23) cr 



3 



<^3, forj^i, 
Ui + P, for j = i, 



(3.24) 



I CTi - ^, /or j = i, 



(3.25) 

P(TOi, TOi+i, CTj, o-,:+i) = e'^*+'''"' - e'^* and 
(3.26) 

(5(r7ii, TOi+i, CTj, o-j+i) = /or f/ie GZRP, 

(3.27) 

P(TOi, TOi+i, (Tj, (T.,+i) = 6°'*+'^"* - e"'* and 
(3.28) 

Q(TOi, m^+i, cr,, cr,:+i) = e-'^'+i - g-'^'+i-'^^'+i for the BLP. 

First notice that ((3l8)) - (l3J9)) and ([3:23)) - ([3:24)) keep ([3lB and ([322]), 
respectively, valid after each jump. Also notice that this theorem reduces to 
Theorem 13.11 in the case of a single shock. In p.25p - p.28p , the exponents 
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CTj + Prrij are simply the parameter values of the Cj marginal or the coupled pair 
(wj, Cj) due to (lOTj) . 

In both the ASEP and the GZRP/BLP, the result can be given a simple 
interpretation. Thinking about the process m{t) of second class particles, the 
right-hand side p.isp is a generator of the interacting system m{t)- Then the 
evolution rules for the densities keep the necessary conditions p.l6p and p.22p 
valid by adjusting the density values to rn{t). 

For the ASEP together with ((X^ - dMU, looks like an asymmetric 

exclusion process except that the (annealed) jump rates of the second class 
particles depend on the density values. Due to the exclusion rule, these particles 
always keep their order, and (j3.20p - (j3.2ip simply say that their right jump rates 
decrease and left jump rates increase as we go from the leftmost to the rightmost 
second class particle. (Recall that p > q and the density increases from left to 
right at each second class particle.) The second class particles thus stay within 
a tight distance from each other, which is already invisible on the macroscopic 
scale of the conservation law. We show in the Appendix that the velocity of 
this structure of shocks agrees with the Rankine Hugoniot velocity (|2.13p taken 
with the leftmost and the rightmost density values. Once the particles, put in 
the shock positions in [5], are replaced with our second class particles, our result 
coincides with the interaction of the random walking shocks of . 

For the exponential GZRP/BLP, ((3?T5)) with the rates ((3?25]) - ((3?28| becomes 
a zero range- type generator i.e., the (annealed) jump rate of a second class 
particle only depends on the local configuration at its site. However, these 
jump rates again depend, via the parameter values a, on the number of other 
second class particles in front and behind. This dependence is such that the 
more behind, the greater right jump rate and smaller left jump rate for a second 
class particle. Lemma 4.3 from [7] is word for word valid also with second class 
particles: the n particles stay within a tight distance from each other, and their 
center of mass performs a drifted simple random walk on the lattice Z/n with 
right jump rate e"^'"" — 6°'"'=" and left jump rate e~°'"8" — e""^'"*' for BLP and 
zero for GZRP. Here aieft is the leftmost, (Xright is the rightmost value of ai. 
The resulting drift of the center of mass of course coincides with the Rankine- 
Hugoniot velocity for a large shock of densities p(crieft) and p(crright)- 

4 Proofs 

4.1 The proof for a single shock 

Proof of Theorem \3.1\ We assume, without loss of generality, that j — 0. To 
prove Theorem 13.11 we take a finite cylinder function Lp on and see how 
the coupled generator (j2.15p acts on under Eq, the expectation w.r.t. j£q. 
This expectation involves summations for all uji variables. We shall change the 
summation variables in each term where the argument of Lp is modified such 
that, after the change, each term will contain ip{ui, C)- 

Lines (|2.16l) - (|2.19l) describe jumps of the second class particle. For all the 
other lines, the second class particle stays at site 0, and the change of variables 
can be done without much complications. In the term with V'(w~^''^7 C"^'*^) in 
the second line, for example, we will put w^^'*^ and ^"^■'^ as our new summation 
variables. This step will change the arguments of the rates, and it will also bring 
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in factors of fi and ft. 

For lines (|2.16p - (|2.19p . however, there is a comphcation. While Eq gives 
probability one on = + 1, this will not hold after the change of variables. 
As an example, the new variables for the term with (^(w^^'", C) in line (I2.16P 
are w"^'*^ and for which we rather have 

C-i = ^Ii'" + 1 and Co — ^a^'^ ■ 

This is the manifestation of the fact that (|2.16p describes a left jump of the 
second class particle from site to site —1. Notice that this new configuration 
(w"^'", C) is now singular to the measure j^q, but not to y__i. Therefore the 
change of variables will also include a change in the expectation from Eq to 
E_i for that term in line (|2.16p . 

The result of changing variables will be of the form 
(4.1) 

Eo(i¥')(^, C) - Eo{^(^^, C)-^}+E-i{^(^, C)-[B+D]}+^^{^{^, 0-[C+E]}, 

where B comes from the term with C) of line (j2.16p . C comes from the 

term with (^(w, C°'^) of line (j2.17L D comes from the term with V5(w, C°'~^) of 
line (j2.18L E comes from the term with <^(w^'°, C) of line (j2.19p . and A from all 
other terms. The aim is to show that A, B + D, and C + E each do not depend 
on the variables uia ■ ■ ■ ^b- In fact, B + D will give Q, C + E will amount to P, 
while A will make up for the negative terms — P — Q in p.2p . We now compute 
the quantities A - E. 
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We start with A and write it line for line as terms come from (|2.15p 

u>.i + l)/^^(a;i+i - 



A r ^ , 1 1^ /(^^« + - 1) , -.1 
^= 1^ + 1, c^.+i - 1) „et,,A„9(,,. ^ P{UJ^, UJ,+i)\ 



i—a — l 
b 



-p(w_i, Wo) +p(w_i, Wo + 1) 

+ p(wo + l,wi-l).^^-^-^^-^ p(wo,wO 



^(wo + 1, wi) +p(wo, wi) 



2^ [g(w, - 1, w,+i + 1) „e(,,A„e(,,. ^ ^(w., ^z+i)J 



2— a— 1 
6 



+ ^[g(w.-l,w.+i + l). ^.(^^)^.(^^^^) 9(w„w.+i)] 

, / 1 , IN - l)A(wo + 1) 

+ q(w_i - 1, Wo + 1) g- — T 9(w-i, Wo) 

M (^^-i)Ai(^o) 

- g(w_i, Wo + 1) + g(w_i, Wo) 

- q(wo, wi) + q{uJo + 1, wi). 

Then 

B + D= [p(w_, + 1, wo - 1) - + 1, wo)] . /(--i + l);(u^o^l) 

(^42) /^(w_i)/i'^(wo) 

, r , , 1^ ^ Nl /(^-i)A(^o) 

+ [^('^-i, Wo + 1) - g(w_i, Wo) ■ —, ^—7 — r, 

^(w_i)pi'^(wo) 

and 

C + E= [Mwo + 1, wi) - p(wo, wO] • 

, [ ( 1 , IN . , Al'^o - 1)m''('^i + 1) 
+ |_g(wo - 1, wi + 1) - g(wo, wi + 1)J fj,0 (cjq) K^i) ' 

We now simplify the expression of A by plugging in the fraction of /z^'s from 
(j2.6p . In the summations we also make use of 

p(wj + 1, LOi+i - 1) • = ^("^^+1; ^i) and 

/(Wi + 1) 

q{uj, - 1, Wi+i + 1) • — — = q(uJi+i, Wi), 

/(,Wj+l + ij 
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a consequence of (|2.4p . and then 

^^^^ UJi) +p(Wi+l, Wi) ~p{uJi, Wj+l) 

= p{uji+i, uji^i) -p(a;j_i, Wj+i) +q{LOi+i, uji^i) - qiuji-i, LOi+i), 
which is a rewriting of (j2.3l) . This latter produces a telescopic sum. 

A = p(a;b+i, UJa-l) - p(Wa-l, Wb+l) + qiuJb+1, ^a-l) - q{(^a-l, (^b+l) 

, f ,1 N e^ ■ fijup - 1) 
+ l)/i(wo) 

+ p(wo + 1, - 1) — — ^ p(wo + 1, wi) 

(4.5) +p(a;o, wi) Wo) 

+ g(a;-i - 1, wo + 1) 77 — ^ ^(w-i, wo + 1) 

+ gl'^-i, '^o) - q{<^n, t^-i) 

+ g(wo, wi + 1) • — — r-T7 —- ~ q[uji, ujq). 

Checking if random walking shocks with second class particles emerge in a 
model now simplifies to checking the existence of a measure /t with which the 
terms A, B + D, C + E each do not depend on uj-i, luq, wi. (Recall that ip 
depends on configurations of sites a . . .b, hence A is allowed to depend on uja-i 
and LOb+i-) When u]™^^ is finite, one also has to take into account the fact that 
B + D, C + E and A are taken under E_i, Ei and Eo, respectively, which give 
zero weight on a;_i = (^max^ _ ^^^max .^^^ _ ^^max^ respectively. 

Substituting the rates of ASEP, ([23]) and give for (g^J, ([t^)) : 

B + D = q-^^Q, C + E=p-\^ = P, 

in agreement with p.Sp . We obtain, also using (j2.8p . from (j4.5p 

A= {p- q){ojb+i - uja-i) - 1- 

Its expectation w.r.t. Eq can be directly computed since (p does not depend on 
uJb+i nor uJa~i, and Eo is product: 

EoA={p-q)-{\-g)-l = -q- \ - p ■ 

1 — A 1 ^ P 

which finishes the proof for ASEP by (gJ]), ([S3) and ([XS]) . 

The rates p(?/, z) = /(y) = e'^^^^^/^^ and q{y, z) = of the exponential 
GZRP result in 

B + D = = Q, C + E = e'^ -e" =P, 

where we used (|2.6p . p. 71) . (|2.9p together with p.6p . This again agrees with 
From (|4.5p we get 

A = /(Wfc+i) - f{uJa-l), 
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which becomes e'^ — e under the expectation of the product measure z/g, also 
in agreement with p.Sp . Similar computations for the exponential BLP lead to 

B + £> = c"'" - c"^ = Q, C + S = c" - c'" = P, 

A = f{0Jb+l) + /(-Wa-l) - f{0Ja-l) - /(-Wb+l), 

which becomes e'' + e~^ — e^ — 6^°" under Eq. □ 
4.2 The proof for multiple shocks 

Proof of Theorem\3M Expectation w.r.t. vf-^^ (|XT4)) will be denoted by W-'^. 
Take a finite cylinder function that depends on pairs (cj^, C,i) ioi a < i < h 
and assume that mi — for i < a and i > b. (This implies, in particular, a 
finite number of second class particles in the system.) We will take the E-— 
expectation of p.9p and change variables in some terms so that after the change 
we have f{uj_, C) in each one. This expectation involves summations for all uji 
variables. We shall change the summation variables in each term where the 
argument of if is modified such that, after the change, each term will contain 
'Pi^ily 0- Lines (|3.10p - p.l3p describe jumps of second class particles. For all the 
other lines, the second class particle stays at site i, and the change of variables 
can be done without much complications. In the term with (^*''+-'^) in 

the first line, for example, we will put and as our new summation 

variables. This step will change the arguments of the rates, and it will also bring 
in factors of and . We will also make use of Q = + mi which 

is a probability one event under E-— . 

For lines (I3.10p - (I3.13p . however, notice that the change of variables influ- 
ences the number of second class particles, and the measures v-^— are singular 
to each other with different vectors m. Therefore with the change of variables 
we will also include a replacement of measures accordingly. 

These changes result in 

b ( 

(4.6) E^'^(i¥.)(c., C) - E C) • 

i=a-\ ^ 

+ E2*^'"-^'^'''{^(^^, C) • [B. + A]} + E2""+^2L-'+^{^(^, Q . + E,\\\^, 

where Bi comes from the term with (/^(a;*''^"'^, C) of line (I3.10p . Ci comes from the 
term with (/5(w, C'^*+^) of line (13. lip . Di comes from the term with (^(w, C*+^^*) 
of line p.l2p . Ei comes from the term with 'p(w*+^'*, C) of line p. 131) . and Ai 
from all other terms. The aim is again to show that Bi + Di, Ci + Ei, and the 
sum of Ai's each do not depend on uja ■■ - UJb. In fact, they will be identified as 
the rates and ((X^ - (jX^ . respectively. We now compute the 
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(4.7) 



quantities Ai - Ei. 

A, = p{uj, + 1, UJ,+i + m,+i - 1 -a m-^ A-a-^im-^i^ \ 

-p{uii, Ui+i + rui+i) 

+ TOi, Wi+i + rrii+i) +p{uj.i, uji+i + m^+i) 
+ glw; + rn^ — 1, uji+i + 1 ■ : ; — ; ^ 

- q{u!i + mi, uji+i) 

- q{uji + rrii, LOi+i + rrii+i) + ^(wi + rrii, uji+i) 

- q{uji, iOi+i) + q{uji + mi, iOi+i) 

= p CJi + 1, W.+i + m,+i - 1 m-^ ^ ^ 

Wi+i) +p{ujt, uj^+i + mj+i) + m^, Wj+i + mi+i) 

+ oIWi +771^ — 1, CJi i 1 + 1 • ; — -, ^ 

- q(uji + TOi, Wi+i + mi+i) + q{uji + mi, cJi+i) - 5(0;^, w^+i), 
+ A = + 1, - 1) - + 1, uji+i + m^+i - 1)] 

+ [q{ui + mi, Wi+i + mi+i) - q{uji + mi, t^i+i)] 

Q + = + mi, w^+i + mi+i) -_p(wi, Wi+i + mi+i)] 

+ [g(a;i - 1, w^+i + 1) - ^(wi + - 1, u^+i + 1)] 



(4.8) 



(4.9) 



X 



+1 



While these formulas are much more general and complicated than (j4.5p . (|4.2p 
and (|4.3p . they can be used to verify Theorem 13.31 
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We start with plugging in the rates of ASEP and (PTTC]) . ([5T7)) . 



Ai = = 0, uji+i = I, rui = 0, mj+i = 0} 



(1 - 



-pl{uji 


= 1, 




hi 


= 0} = 1, uji+1 


= 0, mi+i = 0} 


-pl{uji 


= 0, 




hi 


= 0, nii = 1, TOi+i = 0} 




-pl{uji 


= 1, 




hi 


0, mj = 0, TOi+i = 0} 




+ qi^i 


-1, 




-1 


= 0, mi = 0, mi+i = 0} 


(1 - 




= 0, 




-1 


= 1, mi 0, mi+i = 0} 




- ql{uJi 


= 0, 




-1 


= 0, mi = 0, mi+i = 1} 






= 0, 




-1 


= 1, nii = 0} - ql{uji = 


0, uji+i = 1} 



= {p- q) [a;j+i(l - mj+i) - Wi(l - mi)] 
-pmi(l - mi+i) - (7(1 - mj)mj+i. 

In the second equahty we used that /i^'^(l) = in E-— , hence any situation 
with ujj = 1, mj = 1 does not occur. Validity of the equation can be directly 
checked for the remaining 9 cases of lo' s and m's being or 1. 

Bi + Di^ pl{uji = 0, uji+i = 1, mi = 0, mj+i = 1}- 



+ ^ 0, Wi+i = 0, mi = 0, mi+i = 1} 



1 - g» 

l-g^ 



+1 



= p(l - mi)fn.,+i ■ = (1 - mj)mj+i ■ [(1 - ft)(7 + Qip] 

Qi+l 

using p.l6p and the fact that the indicators imply m*^^'* ~ 1 and therefore 
necessarily Ui = under E^'^^'^'ia'^^ '. Similarly, 

Ci + Ei = pl{ui = 0, Wi+i = 0, mi = 1, mi+i = 0}^^ 

1 - Bi-i 

+ ql{uji = 1, Wi+i = 0, mi = 1, mj+i 0}-^^ 

Qi-l 

= qmi{l - m,j+i) • = mi(l - mi+i) • [(1 - + Qi+iq] . 

Compare these with p.2ip and p.20p . Now we redistribute some terms from 
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the Ai^s and write (|4.6p as 

b 

= 5Z {(1 - - + ftp] 

i=a— 1 

+ TOi(l - TOi+i)[(l - + ft+ig] 

b 

+ (p - g)E2'2i(p(w, C) ^ [uj,+i{l - m,+i) - - m,) 

i—a—l 

- ft+i77ii(l - TTii+i) + gi{l - mi)mi+i]. 

The proof for ASEP is done as soon as we see that the last two lines sum up to 
zero. Without loss of generality, let rrii = for each i < a and i > b. The terms 
with w's are telescopic, and the terms in the bracket will only have uja-i and 
oJb+i in them. These become Qa-i and Qb+i, respectively, under the product 
expectation. Therefore the expectation transforms the summation into 

(4.10) Qb+i-Qa-i- ^ [ft+imi(l - m^+i) - ft(l - mi)mi+i]. 

i—a — l 

Since there are finitely many second class particles, there is a leftmost site ri 
such that = and m^i+i — 1. We successively define 

4 : = inf{i : i > r^, mi ^ 1, m^+i = 0}, and 
Tfc+i : = inf{i : i > r^, mi ^ 0, m^+i = 1}. 

We let inf oo, and K — maxjfc : < oo}. For any i with < i < rk+i we 
have mi — rrii+i — 0, and for any i with rk < i < ik we have to^ = rrii-^-i = 1. 
Therefore, the summand is zero for all these cases, and 

b 

[ft+imi(l - mi+i) - ft(l - mi)mi+i] 

i—a — l 

K K-1 
k=l k=l 

Notice also that by p.l6p Qi is unchanged over an interval without second 
class particles. Hence the sum in the last display is zero, and gii^+i ~ Pb+i, 
Qrx — Qa-1 implies that (|4.10p is zero which finishes the proof for ASEP. 

We now turn to the exponential GZRP and BLP models. Substituting into 
dHZl), gH) and gH) gives 

B, + D, = 0, 
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for the GZRP, and 



,/J((^i+i+m.+i)-/3/2 _ p/J(t^,+m,)-/3/2^ _ ( ^-l3u>,+ i-P/2 _ g-/3a;,-/3/2 



for the BLP, where we also used p.9p . Thus we identified these terms with 
p.25p - (I3.28P as required. Notice that Ai is of gradient form, thus summing it 
in the first line of (|4.6p gives 

b 

i—a— 1 

= /(wf,+i + TOb+i) - f{uja-i + ma-i) for GZRP, 

= fiuJb+1 + TOb+i) - fiuJa-1 + ma-i) - /(-Wfe+i) + /(-Wq_i) for BLP. 

We again can assume finitely many second class particles inside the discrete 
interval [a, b]. Product structure of the measure then allows us to compute the 
above and show that the first line of (14.61) equals 

b 

i—a— 1 

= [e'^o+i - e'""-i] • E^'^(^(cj, C) for GZRP, 

= [e'^'+i - e'^-i - e-'"'+^ + e-"-^] ■ E^'^(^(w, C) for BLP. 

These brackets can be broken up into (notice that ai is constant near the bound- 
aries) 
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E [e---e-]= E [e--e--]= E [' 



i—a—1 i—a—1 

for the GZRP, and 



i+pnii _g-cri + i _|_g-CTi-|-i-^mi-|-i"| 



J2 [e^'+'-e^'-e-^^+'+e-"^] = K'-e"' 

i—a—1 i—a—1 

for the BLP. These give the negative terms in p.lSp and the proof is done. □ 

5 Branching coalescing random walk 
5.1 The model 

In this section we investigate the random walking shock structure found in the 
branching coalescing random walk model (BGRW) by Krebs, Jafarpour, Schiitz 
[5]. This model does not conserve the particle number, it is not a member of 
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the family of Section [5] However, the result we show below is very similar to 
Theorem 13.11 

Besides (|2.1I) . we introduce u/"^ and w'^ corresponding to particle creation 
and annihilation on site i by 

it ) ^fe, for k^i, ji j Wfc, for k ^ i, 

" + for fc = i, " I^Wfe-l, for fc = i. 

The BCRW is a particle system on {0, 1}^, and the dynamics consists of three 
types of processes. Asymmetric nearest neighbour jumps (like in the ASEP): 

cj — > o;*'*"''"'^ with rate p ■ a;i(l — ti^i+i), 

a; — > a;*^"""'* with rate q ■ [1 — 

coalescence (i.e., merging of two particles into one) from the left and from the 
right: 

w — > Lu^^ with rate Cj. • WjWi+i, 

— tj*''" with rate q • WiWi+i, 

and branching (i.e., creation of a new particle next to an existing one) to the 
left and to the right: 

w (jj^"^ with rate bi ■ [1 — 

Lo a/'^^'^ with rate br ■ uji{l — w^+i). 

The factors in the rates here are positive real numbers. The generator of the 
process (with a finite cylinder function ip) is 

(5.1) {L^){ui) = 

b 

^ {p ■ U!i{l - UJ^+l)[ip{ui^-'-^^) - Lp{lil)] + q- {1- UJ.,)LL!i+l[ip{u/^'^''') - (fiiui}] 
i—a~-l 

+ Cr ■ u;iu;i+i[ip{u/^) - ip{uj)] + ci ■ uJiUJi+i[(p{uj^^'^^) - ifiico)] 

+ 6; • (1 - UJi)uJi+l[ip{uj_'^^) - ip{ui)] + br • Wi(l - UJi+l)[ip{ui'+^^) - ifiiuj)]}. 

We will abbreviate B = bi + b,. and C = q + c^. Define, for g G [0, 1], the 
one-site marginal fi^ by 

It is known (and will be apparent from our computations as well) that the 
Bernoulli product distribution 

(5.2) /i^ with the specific density value g* — — — — 

is a translation-invariant stationary distribution. Besides this, there is a trivial 
stationary measure, which is the totally empty lattice with probability one: 
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5.2 Shock measure 



In [5] the following shock measure was considered (adapted to infinite volume): 

(5.3) (g)Ai'^'(g)/. 

This measure is the gi*-stationary distribution on the left of j, and forces the 
empty configuration on the right of j. It has been shown that this shock struc- 
ture performs a biased simple random walk when the algebraic relations (|5.5p 
below hold for the rates. 

This result is similar to that for the ASEP [5] and for the BLP [7j. An 
important difference is however, that whereas in the ASEP and the BLP a good 
definition for the microscopic position of the shock is the position of the second 
class particle, in the BCRW this is not possible (the basic coupling does not 
conserve the number of second class particles). Here, an obvious definition of 
the shock position is the position of the rightmost particle. Indeed, using the 
result of [S], it can easily be shown that the stationary distribution seen from 
the rightmost particle is a product measure with density g* on negative sites 
and density zero on positive ones, which is a similar result to that of [5] and [5]. 
This naturally raises the question whether the rightmost particle itself performs 
a simple random walk in the BCRW. 

Following the previous strategy we define the shock measure 

=(8)M^*(g)Ai'(g)/, 

i<3 i=j i>j 

which is similar to the previous shock measure except that at site j we force 
the presence of a particle. This one is of course the rightmost particle in the 
system. 

Theorem 5.1. The identity 

(5.4) 4:^^■S(^) [a^-^i -m ] + Q- [m. 1 -Ai ] 

^ ' dt—3 t=o — j+i —j' — j-i —r 

holds if 

(5.5) p = 

In this case 

C B C B 

(5.6) P-P-^^ Q = q.^— + c,--^, 

which shows that the rightmost particle performs a simple random walk. 

Notice that (|5.5p does not depend on q, the jump rate of particles away from 
the zero-density region. We remark that both the result in ^ and this theorem 
has a mirror-symmetric version with the empty configuration on the left and 
density g* on the right, p, Cr, and br interchanged with q, ci, and bi, respectively. 

Notice also that 

k— — OG 

exactly equals (j5.3p (just condition on where the rightmost particle is according 
to (|5.3p '). Due to linearity of both sides of (15.41) . our result implies the same 
equation for the measure (|5.3p . thus the result of j8j is recovered. 



B' 
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5.3 Proof for BCRW 

Proof of Theorem \5.1\ Again we set j — 0. The idea is now familiar: first take 
a cylinder function of values uja ■ ■ - ^^b with a<— 1<1<6 and apply the 
expectation Eq w.r.t. //^ of (|5.1|) . Making use of = 1 and Wi = for i > 0, 
this leads to 

-2 

Eo(L(^)(ai) = Eo{ {p■oJ^{l-io,+,Mu'''+^)-^{ul)] 

i—a — 1 

+ c,. • uJiUJi+i[(p{u/^) - (p{uj_)] 
+ ci ■ cj,a;j+i[<y9(a;''+^-'-) - (^(w)] 
+ 5/ • (1 - uji)uji+i[ip{iJ'^) - ip{uj)] 
+ br • a;,(l - - Lp{ui)]] 

+ Cr ■ W_i[(y5(w"^^) - (^(w)] 

+ 6, • (1 -a;_i)[(y5(ti;-i^) -(^(a;)] 

Then change variables to restore ip{uj) in each term, and watch for singularity 
of Ej's for different configurations. In fact and will be taken under Ei, 
0;°'^^ and lJ^^ will be taken under E_i after the change of variables. We arrive 
to 

-2 

(5.7) Eo(L^)(w) =Eo{(^M[ ^ A,+A]}+E_i{(^(w)-i?}+Ei{(^(^)-£;}, 

i—a — l 

where 

Ai^p - {1- uji)ui+i - p ■ ujiil - + q ■ uji{l - cjj+i) - g • (1 - U!i)uji+i 

+ bi ■ ujiUJt+i ■ - 6; • (1 - uji)uji+i 

pe (1) 

+ br ■ WiWi+i ■ — TT-T - Or • W^l - W^+i), 

(1) 
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A ^ -q ■ (l - + Cr • (1 - UJ-l) ■ — TTTTT - Cr ■ UJ-1 - Cl ■ UJ-i 

fie (0) 

M*^' (0) 

+ bl ■ UJ-i ... - 6; ■ (1 - - p-br, 

^le (1) 

D = q he/ — , 

E ^ p- {1 - ujn)—^!--- +br ■ ujQ- 



In the above we freely added or removed the a.s. one factors ujj and 1 — Wj+i for 
terms under Ej. The aim is again to show that D and E become independent 
of the variables and equal to Q and P, respectively, while the A terms amount 
to — Q — P of (j5.6p . Plugging in the specific form (j5.2p of g* simplifies things to 

\p-<i^Cr- — ~bl\ ■ (Wi+i - UJi), 

EC B 
A^ [q- Cr - — -C + bi- — + bi] ■ w_i - q + Cr ■ — - bl - p - br, 



B + C ' P + C" 

B C B 

E = [-p-p - -^+br + br - —] - Wo+p+p - —. 

Since Ai is of gradient form, we can compute the sum in (|5.7I) . We also use 
the fact that ip does not depend on Wa-i, hence this latter has expectation 
g* = B/[B + C]. Then 1^ reads 

C 

Eo(P(y3)(w) = Eo(^(a;){[p- 6^ • 

"^•[pTc + ^]''^-pTc+ p + c '^-^ 

B C B 

+ F,iLp{0l){[-p-p- ■^ + br + br-—] -Uq+P+P- —}. 

Substituting (j5.5p cancels the w-dependence and results in 

C B B + C 

Eo(P<^)(w) =Eo</7(a;){-g- -^-^ -Q • -^-^ ^ P ' ~C~^ 

□ 
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Appendix 



A Shock velocity for multiple shocks in the 
ASEP 

Consider a multiple n-shock in the ASEP. We use a slightly modified notation: 
enumerate the second class particles k — 1 . . .n in spatial order, and denote by 
the constant density value from the position of the k^^ second class particle 
to the left neighbor of the site of the k + second class particle (to infinity for 
k — n). We also denote by g(o) the density value on the left of the I''' second 
class particle. This notation simplifies the arguments to come, and it is possible 
due to the specific fact that second class particles exclude each other and keep 
their order in ASEP. Condition (|3.16l) then becomes 

g(/c)(l - g(fc-l)) _ P_ 

Q(k-i){^ - Q(k)) q 

in this notation, and jump rates (j3.20p . (I3.2ip turn into 
(A.1) 

P(k) = — ■ P = ■ <1 and Q(fc) = — ^ ■ q = ^ p 

1 - ^(fc-i) Q(k~i) 1 - Q(k) Q{k) 

for the fc*'' second class particle. 

First notice that the dynamics of a multiple n-shock, which has been shown 
to be an n-particle simple exclusion dynamics with hopping rates P^j.) (right) 
and Q(k) (left) for particle fc, can be mapped to an n — 1 site open zero range 
process (ZRP), see Figured! The inter-particle distances in the exclusion process 
are interpreted as occupation numbers of the corresponding open ZRP. In this 
ZRP the left and right hop rates from site fc (fc = 1, 2, • • • , n — 1) are independent 
of the occupation number and are equal to Q(k+i) and P(/;) respectively. The 
rate of particle injection on the first and last {n — 1^') site of the system is Q{i) 
and P(„) . 

Due to the attraction of "micro-shocks" in the ASEP, this ZRP has a well- 
defined stationary measure, which has a product structure with site-dependent 
fugacities P2 • The mean velocity of the shock in the ASEP maps to the station- 
ary current in the ZRP. In [35] the stationary current is calculated for open zero 
range processes with constant left and right hopping rates in the bulk. Their 
formula (eq. 15 in that paper) for the current J can easily be generalized for 
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Q(i) P{i) 



Q 



(2) 



P{3) Q(A) P(i) Q{5) P{t 



(5) 
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P(l) Q{2) P{4) Q{5) 

Q{1) P{3) Q{4) P{5] 



o 
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I I I I I 



Figure 2: Mapping of a simple exclusion process (above) to ZRP (below) for 
n = 5. 



inliomogeneous cases as 
(A.2) 



J = 



n Pik) - n Qik) 

71 k — l n 

E n Pn) n Qn) 

k=l 1=1 e=k+i 



Using (jA.l[) we write the denominator of (jA.2|) as 
(A.3) 



k-l 



n h 



1 - 0{e-i) 



k=l i=l ^) l=k+l 

Define now 



9(1.) 6'(o)(l - ^(n)) ^ 



^ = X! £'(fe-i)(l ^ ^ = ^ Q(k-i))Q(k), 

k=l k=l 

for which we have 

p-A = q-B and B - A = £i(o). 
This can be solved to 



and hence (|A.3I) equals 



Q(n) - g(0) ^ 

p-q 

Q{n) - Q{o) g" 



Q(S)){^ - Q{n)) p-q' 

Inserting this into (|A.2p one obtains 



(A.4) J 



em 



g(")~g(o) . 
e(o)(i-e(n)) p-Q 



g(K)(l - g(n)) - g(0)(l - g(0)) 
g(n) ~ g(0) 



One can see that in the rhs. of (|A.4p the usual Rankine-Hugoniot velocity is 
recovered as expected. 
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